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field is calculated in the one-loop approximation of the 2+1 dimensional quantum electrodynamics 
(QED 2 + 1 ) with a nonzero fermion density. We construct the Green function of the Dirac equation 
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I. INTRODUCTION 


Planar charged fermions governed by Dirac equation with external electromagnetic fields attract con¬ 
siderable interest in connection with problems of the quantum Hall effect lj, high-temperature super¬ 
conductivity Q as well as graphene (see, e.g., [1] and 0-0). In graphene, the electron dynamics at low 
energies is described by the two-dimensional Dirac equation for massless fermions 00! though the case 
of massive charged fermions is also of interest Ha¬ 
lt is well known [ll| that every energy level of planar electron in an uniform magnetic field is degenerate 
and the number of degenerate states per unit area is \eB\/(2-Khc) . So, the kinetic energy of the electrons 
can completely quenched at strong magnetic fields. Moreover, the kinetic energy per the Dirac electron 
is of order e « VF\j2\eB\hc (where vf is the dimensionless Fermi-Dirac velocity) and is comparable with 
the Coulomb energy per electron Eq = e 2 /(eo^s), where Ib = \fhcJ\eB\ is the so-called magnetic length 
and eo is the dielectric constant of the medium |6|. In the two-dimensional electron gas, this enhancing 
of Coulomb interactions between electrons in the presence of strong magnetic fields, probably, leads to 
the fractional quantum Hall effect EMI. 

Important physical quantities related to the vacuum polarization are the vacuum charge and current 
densities induced by the background field. Polarization effects in the massive QED 2+1 with a constant 
uniform magnetic field and with a nonzero fermion density were studied in EJE3- In particular, the 
contribution of the induced Clrern-Simons term to the polarization tensor and the effective Lagrangian 
with the electron density corresponding to the occupation of n Landau levels in an uniform magnetic field 
were calculated in fl5| . 

Since the effective fine structure constant in graphene is large, the QED 2+1 effects can be significant 
already in the one-loop approximation. Important quantum relativistic effects were discussed in [16j | (Klein 
paradox) and E3 (Casimir effect). The polarization operator in a strong magnetic held perpendicular 
to the graphene membrane has been calculated in E3| • The problem of light absorption in graphene was 
investigated in fl9j and the Faraday effect in a monolayer graphene sample in a strong constant uniform 
magnetic held perpendicular to it was considered in [20j. 

The induced vacuum current in the held of a solenoid perpendicular to the graphene sample was 
investigated in plj |. and vacuum polarization in the QED 2+1 with an Aharonov-Bohrn (AB) potential 
for massive and massless fermions was studied in [22| . The vacuum electric current due to to vacuum 
polarization in the AB potential for massive case was observed in [23] in “a quantum-tunneling system 
using two-dimensional ionic structures in a linear Paul trap”. Very important phenomenon - charged 
impurity screening in graphene due to the vacuum polarization by a Coulomb held - was investigated in 
[9j, l24H27l | . Effect of spin on the dynamics of the two-dimensional Dirac oscillator in the magnetic cosmic 
string background was considered in [28]. 

In this work, we have calculated the polarization tensor of planar charged fermions in the presence of 
an external constant uniform magnetic held in the one-loop approximation of the QED 2+1 at the hnite 
chemical potential. We have shown that one-loop polarization tensor induces physical effects, which seem 
to be likely to be revealed in a monolayer graphene sample in a strong constant uniform magnetic held 
aligned perpendicularly to the sample. 

We shall adopt the units where c = K = 1. 


II. VACUUM POLARIZATION BY A CONSTANT UNIFORM MAGNETIC FIELD IN 

QED2+1 

The polarization operator (PO) in a constant uniform magnetic held in QED 2+1 is diagonal with respect 
to the photon three-momentum and in the momentum representation is determined by 

/ rl 3 k 

—trWS c (k,B)rS c (k-p,B)}, (1) 

where S c (k, B) is the causal Green function of the Dirac equation with a constant uniform magnetic held 
B in the momentum representation. In the coordinate representation the Green function S c (x fJ/ — x B) 
of the Dirac equation for a fermion of the mass m and charge e in an external constant uniform magnetic 
held in 2+1 dimensions satishes equation 

(- y^Pfj, — m)S(x' J ' — x B) = <5 3 (a; M — x M ), 


( 2 ) 
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where 2 + = x°, x 1 , x 2 =t,x,y is the three-vector, P ^ = — id^ — is the generalized fermion momentum 

operator. The Dirac y^-matrix algebra in 2+1 dimensions is known to be represented in terms of the 
two-dimensional Pauli matrices <jj 

7° = 03 , 7 1 =rrcri, 7 2 = icr 2 , (3) 

where the parameter r = ±1 can label two types of fermions in accordance with the signature of the two- 
dimensional Dirac matrices [29]; it can be applied to characterize two states of the fermion spin (spin ” up” 
and ’’down”) [30j]. We take the magnetic field vector potential in the Cartesian coordinates in the Landau 
gauge A 0 = 0, A\= 0, A 2 = Bx , then the magnetic field is defined as B = diA 2 — d 2 Ai = + 21 , where 

is the electromagnetic field tensor. 

The positive-frequency Dirac equation solutions (the particle states) in the considered field correspond¬ 
ing to the energy eigenvalues (the Landau levels) 


E+ 

= E n = \Jm 2 + 2?z|eB|, n = 0,1,..., 

(4) 

is given by [l5[ 



Vl/ + (t, r) = ^ 

V ^ y/2E n 

( -sign {eE)v£^kL-M ) exp( lE+f + ' nv) • 

(5) 


where the normalized functions U n (z) are expressed through the Hermite polynomials H n (z) as 


U n( Z )= ( 2 n^ !7r l/2)l/2 e ^ ^ H " (*) = VW\(* ~ Vl/eB) 


and p 2 is the eigenvalue because — id y 'S’ + (t, r) = p 2 'b + (t, r). All the energy levels except the lowest level 
(n = 0) with r = 1 for eB > 0 and r = — 1 for eB < 0 are doubly degenerate on spin r = ±1. This 
means that the eigenvalues of the fermion energy except the lowest level are actually spin-independent 
in the configuration under investigation. For definiteness, we consider the case where eB < 0. The 
negative-frequency Dirac equation solutions (the antiparticle states) corresponding to negative energies 
E~ = —E n can be constructed from (J5J by means of the charge-conjugation operation. 

The exact expression for the free electron propagator in an external magnetic field in 3+1 dimensions 
was found for the first time by Schwinger 1311 . We can use the Green function in the considered magnetic 
field obtained for the 2+1 dimensions in [151 ] in the momentum representation in the form 


S C {P) = 


i f dz 

\eB\ J cos z 
0 


exp 


iz 

W\ 


p 0 -m - p“ 


le 


( 7 +r + m) exp(ia 3 z) — 


, tan z 

z 

(7 y+7V) 

cos z 


( 6 ) 


where 2 = \eB\s and s is the “proper time”. The main properties and tensor structure of the PO can be 
obtained from the requirements of relativistic and gauge invariance and also from the symmetry of the 
external field. In the considered external field the PO must be diagonal with respect to the “photon” 
three-momentum and depend only on three independent scalars, which can be constructed from the 
three-momentum £+ and the tensor of external magnetic field F^: 


P 2 = ( V °) 2 - P 2 , Pij.F ijv F vpP p = B 2 p 2 , F^F^ = 2 B 2 . 

We introduce the orthonormalized system of the three-vectors If 

ft = - 4 =( 0 ,-p 2 ) Pi), l 2 = = (p 2 ,PoPhPoP2), I3 = ~^=(po,Pi,P2), 


1 x/p 2 

which satisfy relations 


vVp 2 




= 9ik , £ 9 ik lf If = <T , -{<T - P^/P 2 ) = £ If If, 

i,k j= 1,2 


(7) 

( 8 ) 

(9) 


where g^ is the Minkowski tensor g 11 = g 22 = —g 00 = — 1, and the nonzero diagonal components of g lk 
are g 11 = g 22 = —g 33 = — 1 . 
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The vectors = 1,2 are not eigenvectors of the PO; the PO eigenvalue corresponding to the eigen¬ 
vector is equal to zero due to the gauge invariance 

PlJ n > u ' = w v Pv = 0. (10) 

As a result of calculations, we find the PO in the fully transversal form fl5i ] 


b ) = ^3725172 J —- j dupuw + n 2 IW + C x ier»p p ] x 


c^ 2 dx 


x exp 


ix /pn(l — u 2 ) p 2 sm 1 1 + a /2 smWl - m)/2 2 

— ---:- m + ie 

b \ 4 x sm x 


b = \eB\ 


( 11 ) 


where 


III = Pq (cos ux-u cot x sin ux) — p 2 


w sin x sin 2ua; — 2cosux + coss(l + cos 2 ux) 

• 2 > 

sm x 

II 2 = p 2 (cos ux — u cot x sin ux ), C = 2mcosux (12) 


and e pvp is a fully antisymmetric unit tensor. We note that II pl/ (p) is not symmetric tensor in 2+1 
dimensions. In m last term is the so-called induced Chern-Simons term; using the relation e fMvp p p = 
^ can be written in another form. It should be noted that the mass term in the 
considered QED 2+1 model is not invariant with respect to the operations of spatial (and time) inversion, 
therefore the induced Chern-Simons term must be generated dynamically by the external magnetic field; 
it contributes to the vacuum polarization only in the one-loop QED 2+1 approximation 1~12 . 

In the limit eB = 0, we obtain from OD 




J du[(l — u 2 )(g p, 'p 2 — p p p v ) — 2 ime p ‘ vp p p \ exp 
-1 


P 2 {l-u 2 ) „ 2 , . 

--- to + ie 


13) 


The PO (fldll is a function of only one scalar p 2 and its analytic properties can be studied in the complex 
p 2 plane. It is important that (p) is a real function on the negative real half axis p 2 < 0 what allows 
us to perform integrations in (fl3l) for the domain p 2 < 0 and to obtain [Tsj 


IF» 


e 

47T 




LI V 

p p p 



arctan 





arctan 



(14) 


It should be noted that the free polarization operator in 2+1 dimensions was obtained in another form 
and without the Chern-Simons term in (33]. 

The singularities of IE'' (p) lie on the positive real half axis of p 2 and the point p 2 = 4to 2 is the branch 
point (the threshold for the creation of fermion pairs), so TV ,ll '(p) is an analytic function in the complex 
p 2 plane with a cut [4 to , 00 ); the domains p 2 < 0 and p 2 > 4 m 2 are physical domains, and the domain 
0 < p 2 < 4to 2 is nonphysical. The function II P " (p) in the whole domain of p 2 can be obtained by the 
analytic continuation of (1141) . In the domain p 2 > 4to 2 the free polarization operator gains the imaginary 
part, which is on the upper edge of the cut: 

ImIP» = - py) 4 ”fj-f ■ (15) 

The imaginary part of the PO has a discontinuity 2ImII pl '(p) in going across the cut. 

We can find the polarization operator for charged massless fermions putting in the above formulas 
to = 0. In particular, the free polarization operator for the case to = 0 is a real function on the negative 
real half axis p 2 < 0 and has the extremely simple form: 

e 2 ^U p 2 _ p/ y) 

8 yj— p 2 


rr>,TO = 0) 


(16) 
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The free polarization operator (fTBl) is transverse. Now the point p 2 = 0 is the branch point (the threshold 
for the creation of massless fermion pairs), so II A “ / (p, m = 0) is an analytic function in the complex p 2 
plane with a cut [0, oo). In the domain p 2 > 0 II ^ v {p, m = 0) is pure imaginary and on the upper edge of 
the cut has the form: 


n M1/ (p, to 



p 2 — p^p u ) 

Vp 1 


(17) 


In this form the polarization operator has been calculated in [l9ll33j. In the condensed matter problems, 
the II 00 (p, to = 0) component of polarization tensor is actual that, for example, at p 2 < 0 has the form 


U°°(p,m = 0) 



(18) 


This formula is in a gree ment with that obtained for graphene in the so-called random phase approximation 
in [34| (see, also (aHH). 

It is convenient to represent the polarization operator for the case to = 0 in a weak constant uniform 
magnetic field as follows 


n M "(p, eB, m = 0) = W v {p, m 


0) + e 2 (eB) 2 


tti (p) 



p> l p v \ 


+ TT 2 {p)lili 


(19) 


where the first term is the PO II^"(p, m = 0) and 7Ti (p), 7r 2 (p) are functions of only p. In particular, in 
the domain p 2 < 0 they can be estimated up to constants as 

niip) « C/(-p 2 ) 372 , tt 2 ( p ) « C' lP 2 /(-p 2 ) 5/2 , ( 20 ) 

where C, C\ ~ 0.2. Since fermions are massless dimensionless factors eB/p 2 are built with |p| in place of 

TO. 

To calculate the main contribution in m in a strong magnetic field let us rotate the contour of 
integration in x on — 7r/2 to obtain 


n^(p,b) = --— F e” p ^ 2b (g^ u p 2 — p^p v ) / - —— / G?M(MCotha;sinha; — coshrtx) 

8n 3 / 2 Vb J smh x J 


x exp 


x / pl( 1 - u 2 ) 
b\ 4 


— m 


( 21 ) 


Neglecting the term ~ p\jb in the exponent we integrate (T5T1) in u and obtain 




OO 

, ,,,, 2 U v\ f dx ( COthx 1 \ 


,—xmr/b 


( 22 ) 


We see that the integral is maximum at to = 0, so we can estimate it putting to = 0 in the exponent. But 
the integrand does not contain b at to = 0. Therefore, the leading contribution in b in (1221) is proportional 
to b~ l t 2 


n M "(p, b) 


e 2 C 

4:Tr 3 / 2 y/b 


(g^P 2 


^p v 


C~2. 


(23) 


This result is in agreement with that obtained in [3f>j. Thus, the polarization tensor in QED 2+1 is 
very small (~ |eI3| _1 / 2 ) in a large constant uniform magnetic field. This feature of the PO in QED 2 +i 
essentially differs from that in the massive QED 3 + i, where II 33 component plays a major role and increases 
in a large magnetic field as |37j II 33 ~ |cI3|/to 2 . Physically, this is because only electrons from the lowest 
Landau level couple to the longitudinal components of the photon at \eB\ to 2 [iiclfisl . It is seen that 
the vectors l 1 '', j = 1,2 are spacelike if p 2 > 0, but if p 2 = 0, then the only vector is still spacelike, 

whereas l% —> p^/\fp 2 - It means that the “first” photon mode (~ Z p ) becomes almost free in a very 
strong magnetic field, and the “second” photon mode (~ 1%) therefore does not exist. 
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III. GREEN FUNCTION AND POLARIZATION OPERATOR AT A NONZERO FERMION 

DENSITY 


Now we construct the Green function to the Dirac equation for the case of a nonzero fermion density (the 
finite chemical potential). The Green function with the finite chemical potential can be obtained from the 
Green function in the momentum representation 0 by shifting the variable p° as p° —> p° + p + iSsign(p°), 
where p is the chemical potential. 


S c (p,m,p ) = - 


OO 

i f dz 


exp 


\eB\ J cos z [\eB\ 


iz 


[p° + p + i&ignp 0 ] 2 — to 2 — p 2 ^ + ie 


(7 °{p° + p) + to) exp {ia 3 z) - ^ P + ^ P ^ , S -> +0. (24) 


cos 2 


If p > to, there are real particles occupying the Landau levels, if p < to, then there are no real particles. 
We also assume p > 0 without loss of generality. The integration path passes below the singularities in 
the integrand in (1241) and the imaginary term ifeignp 0 is essential near the poles [39]; the Green function 
(1211) has the poles at the points p° = ±E n — p as well as an imaginary part lmS c (p,p) related to the 
presence of real charged fermions. Rotating the integration path over z in (|24l) into the lower half-plane 
we can extend (1211) onto the whole complex plane of p° with the cuts [p°, 00 ), [— p°_, — 00 ) at the real axis. 
Then, as a function of p° the Green function (12Hl is a limit of some analytic function S c (p°). Denoting 
the integrand in (1211) as S{z), we represent Im S c (p,p) via the discontinuities AS at the edges of the 

OO 

cuts in the form f S{z)dz. For this, we apply method for calculation of the PO discontinuities in the 

— OO 

presence of various external electromagnetic fields in vacuum [37j, l4f| that was extended on the case of a 
nonzero fermion density in [l5j . Finally, we obtain S c (p,p ) in the form 


S c (p,m,p) = -■ 


i f dz 

eB\ J cos z 
0 


exp 


\eB\ 


IZ I r n ■r • m 2 2 2 tan 2 

\p + p + iSsignp] — to — p~ - 


( 7 P + M) - 


—2*7t( 7 P + M), 6 —^ TO. 


(25) 


where 


P° = (po + p) cosz + imsinz, M = tocosz + i(p 0 + p) sin 2 , P = —p/cos 2 , 

P = [to 0 (po + p + m)/2m]5[(po + p - m)/m 0 ]e~ p + 

00 ^ 

—f—777.772.0 Inn (p 2 /|eP|) [<5[(p 0 + d ~ E n )/m 0 } + 5[(p 0 + P + E n )/m 0 }}, 

n= 1 n 

M = [m 0 (po + P + m)/2m]S[(p 0 + p — to) / m 0 ]e _p2 ^ eB '' + 

00 /_1 \n 

+(p° + p)m 0 ^2 P / n „(p 2 /leP|) [5[(p 0 + p - E n )/m 0 \ + <5[(p 0 + P + E n )/m 0 ]} . (26) 

— 1 ^ J ’ n 


Here too is the parameter of the dimension of mass, the Laguerre function I nn (x) = (1 /n\)e~ x ^ 2 L n {x), 
L n (x) is the Laguerre polynomial and all the differences must satisfy inequalities: to + p > 0, E+ — p > 0, 
and —to — p < 0, E~ — p < 0. 

The Green function for charged massless fermions is easily derived from Eqs. dm), m to read 


S c (p,p,m = 0) = 

where 


Ji_ (> + ^ + ? ; ( 5 signp o ] 2 _ p2 ^ + ^ 

—2*7r( 7 P + M), S —y +0. 


i f dz 

- - - / - exp 

\eB\ J cos z 


( 7 P + M) 


(27) 


P° = (po + p)cosz, M = i(po + p)sinz, P = — p/cos 2 , 

P = [to 0 (po + p)/2m]6[(po + p)/m 0 ]\m^oe~ p ' PeB ^ M = [m 0 (p 0 + p)/2m\5[(p 0 +/z)/TO 0 ]|m^oe _p2/|eS| + 
°° ^ ( _ 

+TO 0 (p° + p) Inn (p 2 /|eP|) [<5[(po + P - e n )/m 0 ] + S[(p 0 + p + e„)/m 0 ]], = ±\/2\eB\n. 

£n 


( 28 ) 
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A massless fermion does not have a spin degree of freedom in 2+1 dimensions j4l| but the Dirac equation 
for charged massless fermions in an external magnetic field in 2+1 dimensions keeps the spin parameter. 
Therefore, all the energy levels except the lowest level n = 0 are doubly degenerate; the levels with 
n = n r , t = 1 and n = n r + 1, r = —1 (where n r = 0,1, 2 ... is the radial quantum number) coincide. 

In graphene, the set of Landau levels in an uniform magnetic field aligned perpendicularly to the 
monolayer sample is given by (see, for example, M) 

e n = ±v F ^/2\eB\n, n= 0,1,2,..., (29) 

where v F is the Fermi-Dirac velocity, and the ± signs label the states of positive (electron) and negative 
(hole) energy, respectively. They play the same role as the band index for the conduction (+) and the 
valence (-) band. In addition, the states of positive and negative energy for charged massless fermions 
has the same energy £o = 0 but opposite spins in the ground states. 

We now discuss briefly the polarization tensor (PT) related to contributions coming from real particles. 
We note that the PT contains terms with p ^ 0, B = 0 and p ^ 0, B ^ 0. Terms with p ^ 0, B = 0 are 
very cumbersome and we do not give them. In the absence of magnetic held the n 00 {p, p) component 
of the PT in monolayer graphene has been studied in one-loop approximation in |42 h 44| | . The n 00 (p, p) 
component mainly contribute to the PT in the presence of a weak magnetic held and so it should be 
taken into account in this case. 

Here we give only the expression for the PT related to the combined contribution from real particles 
and magnetic held. It is natural to consider that the medium contains particles, which implies p > m 
and real antiparticles are absent. We also assume that /i ^ m,E n , which means that all the n Landau 
levels are fully filled and no levels are partly filled. 

As a result of long calculations (see, 15|), one can obtain the PT related to the above combined 
contribution in the form 


W v {p,pb) = 


x exp 


27r|m 
(Po + p) 2 


dx[p 2 n r l%l% + mC r x ie plJp pp\ x 


, p 2 (l + * sin 2a: — cos 2a;) 

— 1-—- ex 


2b 


(30) 


where 


_ 1 . 2m(p 0 + v)x nn 
II r — sin + L n 

z b h/ n 

_ i 2m{p 0 + p)x 

Cs r — “ COS-—- + 'iJ-'ri 


2 • 2 
p sin x 

2b 

p 2 sin 2 x 
2b 


sm 


2E n (p 0 + p)x \p 2 - m 2 


cos 


b [ 2\eB\ 

2E n (jp 0 + p)x r p 2 - m 2 
b [ 2\eB\ 


(31) 


Here the hrst terms give the contribution from n = 0 Landau level and [(p 2 — m 2 )/2\eB\\ = N denotes 
the integer part of the function u = ( p 2 — m 2 )/2|eB|, i.e. the largest integer < u. 

The part of the polarization tensor determined by the function n r is of main interest. Calculating it 
for the case of weak constant uniform magnetic held (po + p) 2 |el?|, we obtain 


n Pu (p,pb) = —P 2 n r l p l". 


Here 


n r = 


\eB\ 


(Po + M)[(Po + M) 2 - 4m 2 ] 


6(p- \m\) + 


2|eJ5| 


{Po + ^)[(po + P) 2 - 4 E%\ 


(32) 

xNB(p-E+), (33) 


where d(z) is the Heavyside function; therefore the hrst term gives the contribution if the ground Landau 
level is occupied and the second one contributes if the N Landau levels are occupied. If the magnetic 
held is strong (p 0 + p) 2 \eB\, then only the ground Landau level is occupied in the massive case and 

we obtain 


H r = 


Po + P 
\eB\ ’ 


p > m, |eH| m 2 . 


(34) 


For the case of massless charged fermions we must put m = 0 and replace by in Eq. (l33l) . In a 
strong magnetic held (po + p) 2 /\eB\ <C 1, we also obtain 


n, = ?i±ile W + H%±+xwe( M -4), 


\eB\ 


\eB\ 


(35) 

























where the first and second terms, respectively, give the contributions if the ground Landau level is 
occupied and the N Landau levels are occupied. The main feature is that the polarization tensor (1301) 
with a nonzero fermion density is extremely small (~ \eB\~ 1 11 ) in a strong constant uniform magnetic field. 


IV. DISCUSSION 

We have calculated the polarization tensor in the one-loop approximation of the 2+1 dimensional 
quantum electrodynamics with a nonzero fermion density in a constant uniform magnetic field. The 
polarization tensor contains the contributions from virtual (vacuum) and real charged particles occupying 
the finite number of the Landau levels. In particular, we have found that the polarization tensor in 
QED 2 +i in a strong constant uniform magnetic field is proportional to ~ e 2 \p 2 \/y/\eB\, \eB\ \p 2 \ 
(contribution of virtual particles) and ~ e 2 |p 2 |(po + p)/\sB\,\eB\ (po + p) 2 (contribution of real 
particles). It means that photons become almost free in a very strong magnetic field. 

It should be noted that if one, for instance, needs to investigate the propagation of electromagnetic 
waves in a real graphene strip we must consider it as the three-dimensional object of extremely small but 
nonzero thickness (see, also, 0). 

In addition, we emphasize that the polarization tensors HID and m at nonzero fermion mass are 
finite in the limit p p — > 0 because they contain antisymmetric terms that do not vanish in this limit. The 
coefficient multiplying ie pvp p p in (fill) and (l30l) with p v = 0 is called the induced Chern-Simons coefficient 
and its appearance in the effective Lagrangian of QED 2+1 with an external magnetic field means that 
photons dynamically gain “masses”. The induced Chern-Simons coefficient is calculated exactly and has 
the form 

C C s = —sign(m)|-0(|m| - p) + sign(eB)|- (^0(p - \m\) + N0(p - E+))j , 
where the first term gives the contribution of virtual fermions. 
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